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Molecular abnormalities within the glucocorticoid receptor (GR) stress signaling pathway may confer, or 
reflect, susceptibility to stress in schizophrenia and bipolar disorder, but the extent of such abnormalities in 
the brain is not known. Using RNA-Seq and qPCR in two postmortem cohorts totaling 55 schizophrenia, 34 
bipolar disorder and 55 control individuals, we identified increased FKBP5 and PTGES3 mRNA expression, 
and decreased BAG1 mRNA expression, in the prefrontal cortex in schizophrenia cases relative to controls 
(68.0% [p < 0.001], 26.0% [p < 0.01] and 12.1% [p < 0.05] respectively). We also observed increased FKBP5 
and decreased BAG1 mRNA expression in bipolar disorder (47.5% [p < 0.05] and 14.9% [p < 0.005]). 
There were no diagnostic differences in steady-state FKBP51 protein levels, nor in HSPA1A, HSP90AA1, 
DNAJB1 or HSPB1 mRNA levels. GR, co-factor and chaperone mRNA levels were strongly correlated. 
These results reveal coordinated cortical dysregulation of FKBP5, PTGES3, BAG1 and GR genes within the 
glucocorticoid signaling pathway in psychotic illness. 



Stress has long been considered a contributor to the onset, symptoms and disease course of psychiatric 
illness 18 . In support of this hypothesis, dysregulation of hypothalamic-pituitary-adrenal axis activity, under 
basal conditions and during stress, has been observed in schizophrenia and bipolar disorder 911 . In both 
illnesses, decreased mRNA expression of the glucocorticoid receptor (GR), the primary receptor responsive to 
stress-induced Cortisol secretion, has been observed using multiple post-mortem tissue cohorts and brain 
regions 12-16 . In prefrontal cortex of individuals with schizophrenia and bipolar disorder, increased expression 
of a truncated GR protein isoform has been reproducibly demonstrated 1415 . It is not known whether these 
molecular stress signaling abnormalities in psychotic illness extend beyond GR, to include other cellular com- 
ponents such as stress receptor co-factors and chaperones which are integral to the stress response in the human 
brain. 

Effective GR-mediated stress signaling relies on the coordinated activity of numerous co-factors and chaper- 
ones, which fold GR into a mature state, regulate its affinity for Cortisol, facilitate its translocation into the nucleus, 
and modulate GR-mediated transcriptional activation or repression of target genes (Figure l) 1721 . Some of these 
co-factors and chaperones have been implicated in psychiatric illness. Single nucleotide polymorphisms (SNPs) 
within the gene for FKBP5, which binds GR and facilitates nuclear translocation, have been linked to bipolar 
disorder, major depression and post-traumatic stress disorder 22 26 , and have been associated with impaired 
cognitive performance in individuals with schizophrenia 27 . Increased FKBP5 gene expression has been reported 
in the frontal cortex in individuals with major depression 28 , while decreased FKBP5 mRNA and protein have been 
described in the amygdala in suicide completers 29 . Mood stabilizer treatments for bipolar disorder have been 
shown to up-regulate BAG1 expression 30 . Hsp70 (HSPA1A) and Hsp90 (HSP90B1) gene SNPs have been 
associated with schizophrenia and bipolar disorder respectively 31 " 33 . Increased mRNA expression of Hsp70 
and Hsp27 has been observed in the prefrontal cortex of patients with schizophrenia 34 . Transcriptional changes 
in important stress chaperones suggest their involvement in psychotic illness, yet these findings have yet to be 
replicated, while a number of key molecular regulators of stress responsiveness have yet to be examined. 
Furthermore, since the ability of the brain to respond to stress depends on coordinated modulation of gene 



SCIENTIFIC REPORTS | 3 : 3539 | DOI: 1 0. 1 038/srep03539 



1 



Extracellular space ^.-*" !ow U Cytoplasm 




Figure 1 | The involvement of key co-factors and chaperones in GR-mediated stress signaling. (A) Hsp40 and Hsp70 facilitate, while Bagl impairs, the 
folding of GR into a low steroid-affinity conformation. Bagl also aids degradation of the unstable folded GR complex, (B) Hsp90, p23 and FKBP51 
stablilise the GR complex in a high affinity state, (C) Displacement of FKBP51 by FKBP52 enables the translocation of the cortisol-bound GR 
heterocomplex to the nucleus, (D) Within the nucleus, GR activates or represses the transcription of target genes, at GREs or nGREs respectively. 
Overexpression of Bagl, Hsp90 or p23 can repress GR-responsive gene transcription. Availability of free Cortisol (modulated by CBG) and inter- 
conversion of Cortisol to cortisone by llfj-HSDl/2 also impact GR signaling. See Grad and Picard 17 for comprehensive review. Abbreviations: GRE- 
glucocorticoid response element, nGRE- negative glucocorticoid response element, CBG- Cortisol binding globulin, llfj-HSD- 1 1 (3-hydroxysteroid 
dehydrogenase. 



expression, the relationships between expression of key stress cha- 
perones and co-factors may illuminate how the GR stress pathway is 
regulated in cortical neurons to adapt to the demands of a changing 
environment. 

In this study, we first used RNA sequencing in a cohort of 20 males 
with schizophrenia and matched controls to test if expression of 
stress receptor co-factors and chaperone genes FKBP5, FKBP4, 
PTGES3, BAG1, HSPA1A, HSP90AA1, DNAJB1 and HSPB1 was 
altered in the dorsolateral prefrontal cortex (DLPFC) of people with 
schizophrenia. These genes encode the FKBP51, FKBP52, p23, Bagl, 
Hsp70, Hsp90, Hsp40 and Hsp27 proteins respectively. The DLPFC 
brain region was employed because the most striking diagnostic 
differences in GR mRNA and protein expression in psychiatric ill- 
ness have been evident in that region 1416 . We then used qPCR, in an 
additional independent cohort of up to 35 schizophrenia cases, 34 
bipolar disorder cases and 35 controls, to replicate these results and 
test if any observed changes in stress response-related mRNAs found 
in schizophrenia were shared with, or distinct in, individuals with 
bipolar disorder. Thereafter, informed by our mRNA findings, we 
focused on FKBP5, quantifying FKBP51 protein in schizophrenia 
and bipolar disorder and explored the possible effects of eight 
FKBP5 genetic variants on FKBP5 mRNA and FKBP51 protein in 
the DLPFC. We hypothesized that molecular stress signaling abnor- 
malities in schizophrenia and bipolar disorder may extend beyond 
GR mRNA and protein to include key co-factors and chaperones, 
and that relationships between mRNA transcripts would reflect a 



pathway coordinately dysregulated, resulting in disruption in the 
fine-tuning of stress responsiveness in the human DLPFC. 

Results 

Stress signaling co-factor gene expression in the DLPFC in 
schizophrenia and bipolar disorder. Using RNA-Seq, amplified 
fragments mapping to mRNA transcripts for FKBP5, PTGES3, 
BAG1 and FKBP4 were quantified in the Sydney TRC cohort (for 
cohort details see Table 1). FKBP5 mRNA levels were significantly 
higher in schizophrenia cases than controls (61.5%, ANCOVA F(l, 
33) = 8.93, p = 0.005; Figure 2A). The FKBP5 mRNA expression 
values of 10 (out of 20) schizophrenia cases were higher than the 
highest control FKBP5 expression value. These findings were not 
driven by effects of suicide in the schizophrenia group, because 
although FKBP5 mRNA levels varied significantly according to 
suicide status (ANCOVA F(2, 32) = 6.39, p < 0.005), FKBP5 
expression levels within the schizophrenia (non-suicide) group 
were higher than controls (83% increase), whereas levels in the 
schizophrenia (suicide) group were not different to controls (3% 
decrease). Levels of PTGES3 mRNA were significantly higher in 
schizophrenia cases relative to controls (22.4%, ANCOVA F(l, 33) 
= 7.27, p = 0.01; Figure 2B). There were no significant differences 
between schizophrenia cases and controls in the mRNA expression 
of BAG1 (ANCOVA F(l, 33) = 035, p = 0.56) or FKBP4 (ANCOVA 
F(l, 33) = 1.00, p = 0.33; Figure 2C and 2D). The achieved power of 
ANCOVA analyses ranged from 0.52 to 0.87, and are provided in 
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Age (y) at onset (range) 






23.9(17-40) 


Illness duration (y) ± SD 






23.4 ± 12.41 


Lifetime antipsychotics 






594.94 ± 280.1 1 


(Chlorpromazine equiv., mg) 








Antidepressant use 


yes — 0, no — 20 




yes = 8, no = 1 2 


Type of antidepressant* 






SSRI = 6 , TCA = 2 


Abbreviations: SZ- schizophrenia, M- male, L- left, R- r 


ght, SD- standard deviation, PMI- post-mortem interval, RIN- RNA integrity number, y- years. Data quoted are mean (range) +/— standard deviation. 



Supplementary Table SI. Significant correlations of co-factor 
mRNAs with other clinical/demographic variables in the TRC 
cohort are described in Supplementary Results and Supplementary 
Table SI. 

The number of normalized RNA-Seq reads mapping to specific 
FKBP5 exons in schizophrenia cases and controls were compared 
(Supplementary Results and Supplementary Table S2). Three FKBP5 
exons which are common to FKBP5 mRNA transcripts 1-3, and one 
exon unique to transcript variant 4, had significantly more mapped 
reads in schizophrenia cases than controls (all 1.4 to 1.9-fold 
increases, 0.005 < p < 0.05). 

These results were replicated by qPCR in the Stanley Array cohort 
(for cohort details see Table 2). FKBP5 mRNA levels (total express- 
ion of FKBP5 transcript variants 1-4) differed significantly accord- 
ing to diagnosis (ANCOVA F(2, 87) = 6.00, p < 0.005), with a 68.0% 
increase in schizophrenia cases relative to controls (p < 0.001; 
Figure 2E). As in the TRC cohort, FKBP5 mRNA levels varied sig- 
nificantly according to suicide status (ANCOVA F(4, 85) = 3.67, p < 
0.01), with substantially higher FKBP5 expression levels within the 
schizophrenia (non-suicide) group than controls (87% increase), and 
marginally higher levels in the schizophrenia (suicide) group than 
controls (5% increase). Increased FKBP5 mRNA expression was also 
seen in bipolar disorder cases relative to controls (47.5%, p < 0.05; 
Figure 2E). No effects of suicide were seen in the bipolar disorder 
group. The FKBP5 mRNA expression values of 10 (out of 32) schizo- 
phrenia cases, and 6 (out of 29) bipolar disorder cases, were higher 
than the highest control FKBP5 expression value. PTGES3 mRNA 
expression was also 35.8% greater in schizophrenia cases than con- 
trols (Figure 2F). When testing whether findings from the TRC 
cohort were replicated, this increase was significant (one-tailed t-test, 
p < 0.01), although the overall effect of diagnosis on PTGES3 
expression did not reach significance (ANCOVA F(2, 88) = 2.70, 
p = 0.07). Unlike in the TRC cohort, BAG1 mRNA levels differed 
significantly between diagnostic groups (ANCOVA F(2, 91) = 5.02, 
p < 0.01), with a 12.1% decrease in schizophrenia cases relative to 
controls (p < 0.05), and a 14.9% decrease in bipolar disorder cases 
relative to controls (p < 0.005; Figure 2G). There were no significant 
diagnostic differences in mRNA expression of FKBP4 (ANCOVA 
F(2, 91) = 2.11, p = 0.13; Figure 2H). The achieved power of 
ANCOVA analyses, and significant correlations of co-factor 
mRNAs with other clinical/demographic variables in the Stanley 
Array cohort ranged from 0.74 to 0.83, and are provided in 
Supplementary Table SI. 



Stress signaling chaperone gene expression in the DLPFC in schi- 
zophrenia and bipolar disorder. When stress signaling chaperones 
were quantified by RNA-Seq in the Sydney TRC cohort, there were 
no significant differences between schizophrenia cases and controls 
in their mRNA expression levels: HSPA1A (ANCOVA F(l,33) = 
0.01, p = 0.99), HSP90AA1 (ANCOVA F(l,33) = 2.85, p = 0.10), 
DNAJB1 (ANCOVA F(l,33) = 0.02,p = 0.88) orHSPBl (ANCOVA 
F(l,33) = 0.38, p = 0.54) (Figure 3A-D). Similarly, in the Stanley 
cohort, there were no significant differences in mRNA levels of 
HSPA1A (ANCOVA F(2, 89) = 0.10, p = 0.90), HSP90AA1 
(ANCOVA F(2, 89) = 0.69, p = 0.50), DNAJB1 (ANCOVA F(2, 
85) = 0.11, p = 0.89) or HSPB1 (ANCOVA F(2, 90) = 0.73, p = 
0.48) between schizophrenia, bipolar disorder and control groups 
(Figure 3E-H). 

Relationships between mRNA expression of GR transcripts, co- 
factors and chaperones in the DLPFC in the Stanley Array cohort. 

Spearman's correlations, using non-transformed data, were 
performed to explore the relationships between mRNA expression 
of GR transcripts, co-factors and chaperones in the DLPFC 
(Figure 4). Expression of total GR, GR-1B, GR-1C and GR-1H 
mRNA transcripts was previously quantified in the same cohort 
using the same qPCR methods 16 . The HSP chaperones and 
PTGES3 were significantly positively correlated with each other 
(all 0.29 < r < 0.77, p < 0.01), as were co-factors BAG1 and 
FKBP4 (r = 0.45, p < 5xl0~ 6 ). GR-1B mRNA was negatively 
correlated with FKBP5, PTGES3 and all HSPs (all -0.19 < r < 
-0.53, lxlO" 6 < p < 0.08) and positively correlated with FKBP4 
and BAG1 (all 0.23 < r < 0.33, p < 0.05). FKBP5 was positively 
correlated with PTGES3 and all HSPs (all 0.20 < r < 0.29, p < 0.06), 
and negatively correlated FKBP4 and BAG1 [trend-level 
significance: r = -0.18 (p = 0.09) and r = -0.19 (p = 0.06) 
respectively]. FKBP4 and BAG1 were not significantly correlated 
with any HSPs. 

FKBP51 protein abundance in the DLPFC in schizophrenia and 
bipolar disorder. To determine whether FKBP5 mRNA alterations 
in schizophrenia and bipolar disorder (involving transcript variants 
1-4) were accompanied by changes at the protein level, FKBP51 
protein was quantified in the Stanley Array cohort (Figure 5A), 
using an antibody targeting FKBP51 isoform 1, which is translated 
from mRNA transcript variants 1, 2 and 3. There was no significant 
effect of diagnosis on FKBP51 protein abundance (ANCOVA F(2, 
95) = 0.96, p = 0.39; Figure 5B). 
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Figure 2 | Stress signaling co-factor gene expression in the DLPFC. (A-D) Normalised expression of FKBP5, PTGES3, BAG1 and FKBP4 genes in 
schizophrenia cases (n = 20) and controls (n = 20) from the Sydney TRC cohort; (E-G) Normalised expression of FKBP5, PTGES3, BAG1 and FKBP4 in 
schizophrenia cases (n = 35), bipolar disorder cases (n = 31) and controls (n = 34) from the Stanley Array cohort. * p < 0.05, **p < 0.005. 
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Table 2 Demographics of 


schizophrenia, bipolar disorder and matched control subjects in the Stanley Array post-mortem brain tissue 


cohort used in this study 








Stanley Array cohort 


Control group (n — oO) 


ochizophren ia group (n — Jo) 


Dipolor disorder group (n — 


Diagnostic subtype 


_ 


SZ(disorg.) = 1 SZ(paranoid) = 8 


BPl = 27, BP2 = 4, BPNOS = 2, 






SZ(undiff.) = 26 


schizoaffective = 1 


A 

Age in years (range) 


44.2 (31-60) 


42.6 (19-59) 


45.4 (19-64) 


Gender 


9F, 26M 


9F, 26M 


1 8F, 1 6M 


Hemisphere 


1 6L, 1 9R 


1 71, 1 8R 


19L, 15R 


pn — ou 


6.61 +/- 0.27 


6.48 +/- 0.24 


6.43 +/- 0.30 


PKA\ /k^i irc\ -t- 

r/vti (noursj n ou 


29.4 +/- 12.9 


31.4 +/" 15.5 


37.9 +/- 18.6 


dim + en 

F\l 1 N — JLJ 


7.23 +/- 0.87 


7.36 +/- 0.61 


7.34 +/- 0.88 


Manner of death 


natural = 35 (heart 


natural = 28, suicide = 7, (heart attack 


natural = 1 9, suicide = 1 5 (heart attack 




attack = 31, other = 4) 


= 15, overdose = 4, pneumonia = 5, 


= 7, overdose = 1 2, drowning = 3, 




hanging = 3, jumped = 3, other = 5) 


asphyxiation = 2, hanging = 2, 








other = 8) 


Age (y) at onset (range) 




21 .3 (9-34) 


25.3 (14-48) 


Illness duration (y) ± SD 




21.3 +/- 10.2 


20.2 +/- 9.6 


Lifetime antipsychotics (Fluphenazine 


85004 +/- 100335 


10212 +/-22871 


equiv., mg) 








Antidepressant use 


yes = 0, no = 35 


yes = 9, no = 26 


yes = 1 9, no = 1 5 


Type of antidepressant* 




SSRI = 4SNRI = 0SARI = 2, 


SSRI = 9 SNRI = 4 SARI = 5 TCA = 6, 




TCA = 2, other = 2 


other = 1 


Abbreviations: M- male, L- left, R- right, SD- standard deviation, PMI- post-mortem interval, RIN- RNA integrity number, y- years, BP1- bipolar disor 


der type 1, BP2- bipolar disorder type II, BPNOS- bipolar 


disorder not otherwise specified, SZ- sch 


izophrenia. Data quoted are mean (range) +/- 


standard deviation. *Note- some individuals took multiple antidepressant medications. 



Effects of FKBP5 genotype on FKBP5 mRNA and FKBP51 pro- 
tein expression in the DLPFC. Individuals in the Stanley Array 
cohort were genotyped at 8 SNPs within the FKBP5 gene (Figure 5C). 
A significant main effect of genotype on FKBP5 mRNA levels was 
observed for rs4713916 (ANOVA F(2, 76) = 3.49, p < 0.05) within 
combined control, bipolar disorder and schizophrenia samples. 
FKBP5 mRNA expression was 27.5% lower in TT homozygotes (n 
= 9) than TC heterozygotes (n = 32, p < 0.05). There were no other 
significant main effects of genotyped SNPs on FKBP5 gene 
expression (all F < 2.23, p > 0.11), nor were there significant 
genotype x diagnosis interactions. There were also no significant 
main effects of genotype on FKBP51 protein levels in the DLPFC 
(all F < 1.77, p > 0.18). 

Discussion 

We have identified, for the first time, alterations in FKBP5 transcrip- 
tional regulation in the brain in schizophrenia and bipolar disorder, 
consistent with accumulating evidence for a role of the FKBP5 gene 
in the etiology of psychiatric illness 23 ' 24,26,28 . We also provide evidence 
for the involvement of two other GR co-factors, PTGES3 and BAG1, 
in psychotic illness. Considering these findings in the context of our 
previous work with the Stanley Array schizophrenia/bipolar disorder 
cohort, it is likely that altered expression of FKBP5, PTGES3 and 
BAG1 mRNAs does not occur independently, but occurs in parallel 
with decreased expression of multiple GR mRNA transcripts 16 and 
increased expression of the functional, truncated GRcx-Dl protein 
isoform 35 . Taken together, our results suggest widespread stress sig- 
naling alterations in both schizophrenia and bipolar illness. 

We did not replicate the previous report of increased HSPA1A, 
HSPA1B and HSPB1 expression in the DLPFC found in 14 patients 
with schizophrenia relative to controls 34 . This may be due to the large 
variance in mRNA observed for these transcripts perhaps due to 
increased biological heterogeneity observed in schizophrenia, or dif- 
ferences in the sensitivities and target specificities of microarray, 
RNA-Seq and qPCR methods. Diagnostic differences were restricted 
to co-factors (FKBP5, PTGES3, BAG1) whose functions more spe- 
cifically relate to steroid hormone signalling, and were not seen in 
HSP chaperones which have a diverse range of additional functions. 

FKBP51 regulates stress responsiveness by inhibiting GR nuclear 
translocation (Figure l) 17 . However, given that increases in FKBP51 



protein were not observed in schizophrenia and bipolar disorder, it is 
somewhat unclear what functional impact mRNA alterations may 
have in the DLPFC. In this study, RNA-Seq mRNA alterations arose 
from both FKBP5 mRNA variants 1-3 collectively, and FKBP5 
mRNA variant 4 (see Supplementary Table S2), while qPCR changes 
were detected using a pan FKBP5 mRNA probe. At the protein level, 
we measured FKBP51 isoform 1, which is translated from FKBP5 
mRNA transcripts 1-3, but not the truncated FKBP51 isoform 2, 
which is translated from FKBP5 mRNA transcript 4. Therefore, it 
is possible that we failed to detect diagnostic differences in FKBP51 
isoform 2, which arose from dysregulation of FKBP5 mRNA variant 
4 expression. For FKBP51 isoform 1, since we measured steady-state 
protein levels (one snap-shot in time), and this molecule is involved 
in many active cellular process, it is possible that a large increase in 
FKBP5 mRNA would result in more FKBP51 protein that is then 
readily used and degraded, resulting in no apparent change in overall 
protein levels. Further studies are needed to clarify if FKBP5 turn- 
over, subcellular localization or binding to GR is changed in people 
with psychotic mental illness. A discordance of FKBP5 mRNA and 
protein findings has also previously reported in the context of geno- 
type effects in healthy individuals 23 . The extent of FKBP5 mRNA 
upregulation in psychotic illness observed in this study, and its 
occurrence in two independent post-mortem cohorts, suggest that 
FKBP5 mRNA over-expression may be a robust feature of these 
illnesses. 

It is possible that FKBP5 mRNA over-expression in the DLPFC 
may occur in schizophrenia and bipolar disorder cases due to hyper- 
cortisolemia, which has been described in these illnesses 10,11 . FKBP5 
gene expression is upregulated by glucocorticoids through GR in 
rodent and primate brain 36-40 , and in squirrel monkeys, chronic 
hypercortisolemia is compensated for by over-expression of 
FKBP5 41 . Therefore, it is possible that FKBP5 upregulation may have 
arisen in an attempt to 'buffer' chronically increased stress/cortisol in 
individuals in this cohort. Consistent with this hypothesis, FKBP5 
mRNA was negatively correlated with GR mRNA, which is known to 
be downregulated by hypercortisolemia 42 . Unfortunately, data on 
circulating Cortisol levels for individuals in these cohorts are not 
available. Alternatively, 'high-induction' FKBP5 polymorphisms 
may contribute to FKBP5 mRNA upregulation in psychiatric illness. 
The rsl360780 FKBP5 polymorphism has been shown to influence 
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Figure 3 | mRNA expression of HSP chaperones in the DLPFC. (A-D) Normalised expression of HSPA1A, HSP90AA1, DNAJB1 and HSPB1 genes in 
schizophrenia cases (n = 20) and controls (n = 20) from the Sydney TRC cohort; (E-G) Normalised expression of HSPA1A, HSP90AA1, DNAJB1 and 
HSPB1 in schizophrenia cases (n = 35), bipolar disorder cases (n = 31) and controls (n = 34) from the Stanley Array cohort. 
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FKBP5 gene transcription in vitro 26 , and FKBP5 mRNA/FKBP51 
protein expression in blood lymphocytes 23 . Neither this SNP, nor 
the seven others genotyped, were over-represented in schizophrenia 
and bipolar disorder cases in this study, and therefore are not likely to 
contribute to the observed diagnostic mRNA differences. Albeit with 
small sample size for genetic association, our data do provide some 
supporting evidence for a biological effect of FKBP5 polymorphisms 
in human frontal cortex, since we found that the rs4713916 TT 
genotype was associated with reduced FKBP5 mRNA levels in the 
DLPFC. This finding supports previous indications that the TT geno- 
type may be functional, as it is linked with impaired recovery from 
social stress 43 . How impaired stress recovery may relate mechanisti- 
cally to decreased prefrontal FKBP5 mRNA expression, as observed 
in TT homozygotes in this study, remains to be determined. 

P23, encoded by the gene PTGES3, also plays a key role in gluco- 
corticoid signaling and was increased at the mRNA level in the 
DLPFC in individuals with schizophrenia in our study. In the GR 
heterocomplex in vitro, p23 is an obligatory co-factor 19 and is the 
limiting component of the complex 44 , functioning to stabilize the 
interaction of the complex with GR in the cytoplasm. Paradoxi- 
cally, p23 also acts in the nucleus to inhibit GR-mediated gene tran- 
scription 45 and disassemble GR transcriptional machinery in the 
nucleus 46 . In vivo, p23 is critical for appropriate glucocorticoid 
responsiveness 47 . As a result, any increase in p23 protein synthesis, 
abundance and/or function in schizophrenia, which could occur due 
to increased PTGES3 mRNA expression, could profoundly impact 
cellular stress signaling. 

Downregulation of BAG1 gene expression in schizophrenia and 
bipolar disorder is consistent with other recent work that also impli- 
cates Bagl in the pathophysiology of affective illness. A significant 
decrease of small magnitude (12.1%) was observed in the Stanley 
Array cohort, but not in the TRC cohort, possibly due to the smaller 
TRC cohort sample size, or neurobiological differences between indi- 



viduals in the two cohorts. BAG1 impairs GR folding by interacting 
with Hsp70 4S , facilitates complex degradation 49 and inhibits GR- 
mediated target gene transactivation 50,51 . Upregulation of BAG1 
due to chronic stress has been observed selectively in female rats 52 . 
Mood stabilizer treatments for bipolar disorder have been shown to 
upregulate BAG1 expression 30 , while BAG1 downregulation may 
impair recovery from stress-induced depression-like behaviors 53 . 
Although BAG1 mRNA did not vary according to antidepressant 
or antipsychotic exposure in this study, it is possible that decreased 
BAG1 expression may be linked to affective disturbances that can 
occur in both schizophrenia and bipolar disorder. 

Strong correlations between mRNA expression of GR, co-factors 
and chaperones were evident in this study. Counter- intuitively, the 
mRNA expression of GR and 4/5 of its obligatory co-factors/chaper- 
ones (PTGES3, DNAJB1, HSPA1A and HSP90AA1) were inversely 
correlated. It is plausible that these co-factors and chaperones are co- 
regulated by a common mechanism, whose function is to normalize 
stress signaling by compensating for changes in GR expression, cre- 
ating homeostasis. Interestingly, the expression of two co-factors, 
FKBP5 and FKBP4, which inhibit and facilitate GR nuclear transloca- 
tion respectively, were also negatively correlated, suggesting that these 
co-factors are not coordinatedly up-regulated to 'balance' their effects, 
but may be reciprocally regulated to exaggerate their individual effects. 

In this study, we have identified dysregulation of stress signaling 
cofactors FKBP5, PTGES3 and BAG1 in the DLPFC of bipolar dis- 
order and/or schizophrenia patients, within a network of transcripts 
(consisting of GR, HSP chaperones and cofactors) that were highly 
correlated. We hypothesize that the coordinated regulation of stress 
signaling transcripts is integral to appropriate stress responsiveness 
in the frontal cortex, and that this process may be disrupted in those 
with psychiatric illness. Future work investigating GR, FKBP5, 
PTGES3 and BAG1 in schizophrenia and bipolar disorder may 
identify mechanisms underlying failed adaptation to stress in schizo- 



SCIENTIFIC REPORTS | 3 : 3539 | DOI: 10.1038/srep03539 



7 



A 



75 kDa 
50 kDa 



actin 



STD BP BP C BP C BP SZ C C SZ SZ BP 





Location 


position of gene on 


SNP 








(UCSC build 


chr. 6, relative to 


(major/ 




HWE 


dbSNP 


Hgl9, Feb 


FKBP5 


minor 




P- 


rs# 


2009) 


transcript variant 1 


allele) 


MAF 


value 


rs4713916 


35669983 


promoter region 


C/T 


0.306 


0.388 






(upstream of exon 1 ) 








rs9470080 


35646435 


intron 1 


G/A 


0.350 


0.442 


rs47 13902 


35614026 


intron 1 


A/G 


0.237 


0.430 


rs 1360780 


35607571 


intron 2 


G/A 


0.330 


0.616 


rs9296158 


35567082 


intron 5 


C/T 


0.330 


0.616 


rs7757037 


35548236 


intron 8 


C/T 


0.445 


0.465 


rs3800373 


35542476 


3'UTR 


T/G 


0.305 


0.549 


rs 1043805 


35541432 


3' UTR 


T A 


0.177 


0.949 




CC TC 

rs4713916 genotype 



TT 



Figure 5 | FKBP51 protein expression and FKBP5 genotype effects in the Stanley Array cohort. (A) representative FKBP51 western blot; (B) FKBP51 
protein abundance in schizophrenia, bipolar disorder and control groups; (C) Details of genotyped FKBP5 SNPs. For each SNP, the bases occurring on 
the minus strand are described; (D) FKBP5 mRNA levels in the DLPFC according to rs47 139 16 genotype, in CC homozygotes (n = 44), TC heterozygotes 
(n — 32) and TT homozygotes (n — 9). Abbreviations: C- control, BP- bipolar disorder, SZ- schizophrenia, kDa- kilodaltons, MAF- minor allele 
frequency, HWE- Hardy- Weinberg Equilibrium. 



phrenia and bipolar disorder, and pave the way for strategies to better 
buffer the effects of stress in individuals vulnerable to poor stress 
responses. 

Methods 

Tissue collection. This study was carried out in accordance with the declaration of 
Helsinki, after approval by the Human Research Ethics Committee at the University 
of NSW (#HREC07261). Written informed consent for use of tissues in the study was 
obtained from next of kin. 

RNA sequencing study. Tissue from the DLPFC (BA46) of 20 controls and 20 
individuals with schizophrenia (the Sydney TRC cohort) was obtained from the New 
South Wales Tissue Resource Centre (Table l) 54 . The RNA-Seq performed in this 
study, employing the SOLiD system (Life Technologies, Carlsbad, CA) has been 
previously described in detail 55 . On average, —50% of the 135,355,990 unpaired reads 
per individual, each 50 nucleotides in length, were mapped to the genome, consistent 
with other studies 56 . There were no diagnostic differences in read length or average 
read count. The mapped reads were aggregated into 33,105 transcripts. For each 
individual sample, an aggregate expression value for each co-factor and chaperone 
gene was calculated by taking the geometric mean of the number of mapped reads 
(normalised for total read length and total number of sequencing reads) for all known 
transcripts of that gene. Further information is available as Supplementary Methods. 

Quantitative real-time PCR, western blotting and genotyping study. DLPFC 
tissues (middle frontal gyrus) from the Stanley Array collection cohort were used for 
qPCR, western blotting and genotyping (Table 2). Genomic DNA (from 34 
schizophrenia/30 bipolar disorder/31 control), total mRNA (from 35 schizophrenia/ 
31 bipolar disorder/34 control), and crude protein homogenate (from 35 
schizophrenia/34 bipolar disorder/35 control) were provided by the Stanley Medical 
Research Institute, having been extracted as previously described 15 . 

Pre-designed Taqman Gene Expression Assays (Applied Biosystems, Foster City, CA, 
USA) were chosen for each of the genes analysed: FKBP5 (Hs01561006_ml), FKBP4 
(Hs00427038_gl), PTGES3 (Hs00832847_gH), BAG1 (Hs00185390_ml), HSPAlA 
(Hs00271229_sl), HSP90AA1 (Hs00743767_sH), DNAJB1 (Hs00428680_ml) and 
HSPB1 (Hs00356629_gl). For FKBP5, the Taqman primer/probe set detected all 
FKBP5 mRNA transcript variants (1-4), at the exon 5-6 junction (variants 1, 3 and 4) 



or exon 6-7 junction (variant 2). All reactions were performed in triplicate. Expression 
levels were normalized to the geometric mean of four validated 'housekeeper' genes, 
from distinct biological pathways, whose raw expression levels did not differ between 
diagnostic groups: ACTB (Hs99999903_ml), B2M (Hs99999907_ml), TBP 
(Hs00427620_ml) and UBC (Hs00824723_ml). 

Western blotting was conducted and data normalized as previously described 14 , 
using the ab2901 anti-FKBP51 primary antibody (Abeam, Cambridge, UK), which 
targets amino acids 395-407 of FKBP51 isoform 1 (translated from FKBP5 transcript 
variants 1-3). Samples were run in duplicate in separate experimental runs, and the 
geometric mean of both runs calculated. Further information is available as 
Supplementary Methods. 

Eight promoter or putative functional FKBP5 gene SNPs were genotyped 
(Figure 5) at the Australian Genome Research Facility using Sequenom MassArray 
method previously described 16 , with an assay pass rate of 96-98%. PLINK software 
(version 1.06, http://pngu.mgh.harvard.edu/purcell/plink) 57 was used for Hardy- 
Weinberg equilibrium testing. 

Prior to analysis, all data were either approximately normally distributed, or log 
transformed to yield an approximately normal distribution (skewness between — 1 
and 1). Population outliers, whose values were greater than 2 standard deviations 
from the diagnostic group mean, were excluded from all analyses. For all mRNA 
measures, Pearson correlation analyses were conducted with age, brain pH, PMI and 
RIN for all individuals. Within schizophrenia and bipolar disorder cases, correlation 
analyses were also conducted with age-of-onset, duration-of-illness and antipsychotic 
drug measures (ascertained by converting all antipsychotic drugs taken to lifetime 
chlorpromazine [TRC cohort] or fluphenazine [Stanley Array cohort] equivalents). 
For all mRNA transcripts, analyses of covariance (ANCOVA) were used with age, pH 
and RIN as covariates, to identify effects of diagnosis and gender. For analysis of 
suicide and antidepressant use, schizophrenia and bipolar disorder groups were 
subdivided according to suicide and antidepressant status) and ANCOVA analysis 
performed. For FKBP51 protein, ANCOVA was used with age, pH and PMI as 
covariates. The observed power of all ANCOVA analyses is given in Supplementary 
Table SI. Planned comparisons of schizophrenia, bipolar disorder and control groups 
were made using Fisher's LSD post-hoc test. Two-way ANOVAs were also used, with 
diagnosis and genotype as independent factors, to identify the effects of genotype on 
mRNA expression. To analyze correlations between GR transcripts, co-factors and 
chaperones, Spearman's correlations with non-transformed data were used. Statistica 
software (Statsoft, Oklahoma, USA) was used for all analyses. 
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